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Abstract. For conventional superconductors the electron-electron interaction is 
strongly reduced by retardation effects, making the formation of Cooper pairs pos- 
sible. In the alkali-doped Fullerides, however, there are no strong retardation effects. 
But dielectric screening can reduce the electron-electron interaction sufficiently, if we 
assume that the random-phase approximation (RPA) is valid. It is not clear, however, 
if this assumption holds, since the alkali-doped Fullerides are strongly correlated sys- 
tems close to a Mott transition. To test the validity of the RPA for these systems we 
have calculated the screening of a test charge using quantum Monte Carlo. 



INTRODUCTION 

In order to lead to an effective attraction between electrons, the electron-phonon 
interaction has to overcome the strong electron-electron repulsion. For conventional 
superconductors retardation effects strongly reduce the Coulomb repulsion [1]. The 
resulting effective interaction is described by the dimensionless Coulomb pseudopo- 
tential //*. In the above sense, the alkali-doped Fullerides are non-conventional 
superconductors, since for them retardation effects are very inefficient in reducing 
the Coulomb repulsion [2]. One can, however, conceive an alternative mechanism 
for reducing /!*: Within the random phase approximation (RPA) the Coulomb 
interaction is strongly reduced by the very efficient dielectric screening [2]. 

It is not clear, however, how well the random phase approximation describes 
the screening in the alkali-doped Fullerides. RPA works well in the limit of weak 
interaction, while in the opposite limit it is qualitatively wrong. Very little is known 
about the screening in the intermediate region, where interaction and kinetic energy 
are comparable. This is the region where the superconducting Fullerides, being 
close to a Mott transition, lie. 

To find out how well RPA works for the alkali-doped Fullerides, we have stud- 
ied the static screening of a test charge in a Hubbard-like model using quantum 
Monte Carlo methods. Our results indicate that RPA gives a surprisingly accurate 
description of the static screening on the metallic side of a Mott transition. 



QUANTUM MONTE CARLO CALCULATIONS 



To describe the conduction electrons in A 3 C 60 we use the multi-band Hubbard- 
like Hamiltonian 

Hq ^ ' tij mm 1 C^ rna Cj m i a -\- U ^ t f^ima^im' a' ■ 

m,m',a i, mo-Cm' a' 

The indices i, j label the sites of the fee lattice, denoting nearest neighbors. 
Each Cgo-molecule randomly takes one of two orientations [3] . The hopping matrix 
elements tij mm ' between the three-fold degenerate t\ u orbitals on neighboring sites 
are obtained from a tight-binding parameterization [4]. The Hubbard interaction U 
is varied to study the effect of correlations. Experimental estimates give U ~ 1.5 eV 
[5]. 

To study the static screening, we introduce a test charge q on a site i q . The 
corresponding term in the Hamiltonian is 

Hi(q) = qUj2 n i g ma- 

ma 

We can then find the response of the system upon introduction of the test charge 
by comparing the electron density n q on site i q for the system with the test charge, 
to the electron density no for the unperturbed system: An = no — n q . 

We determine the ground state properties of the Hamiltonians H and H + Hi 
(for finite fee clusters) using diffusion Monte Carlo with the fixed-node approx- 
imation [6]. In this method we construct a trial wave function \^t) and allow 
it to diffuse towards the exact (within the fixed node approximation) solution 
IV'o)- F° r the trial function, we make a generalized Gutzwiller ansatz: (R\ty T ) = 
gd(R)gn q (R) |<|>} i s the Slater determinant calculated using the random phase 

approximation. The Gutzwiller factors reflect the two interaction terms, where 
for a given configuration R of the electrons d(R) = (R\J2i,mo<m'o> n imanim'a'\R) 
and n q (R) = (R\ J2ma n i q ma\R) ■ We optimize the free parameters g and g q using a 
correlated sampling technique. 

Diffusion Monte Carlo only provides us with the mixed estimator 
(^t\0\^o) / (^t\^o) °f an observable O. If, however, the trial function |^ T ) is 
close to the ground state, we can estimate the ground state expectation value using 
the extrapolated estimator 

(^o|C|^o) _ 2 (#t|£>|*o) _ (^t\0\^ t ) 

The last term, the expectation value of the observable for the trial function, is easily 
calculated using variational Monte Carlo. To make sure that the extrapolated esti- 
mator gives reliable results, we have compared our Monte Carlo calculations with 
the results from the exact diagonalization for small systems. For larger systems, 
where exact diagonalization is not possible, we have checked that the trial function 
is close to the ground state wave function. 
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FIGURE 1. Screening of a test charge q = 0.5 e as a function of U for a cluster of 32 molecules 
(band width W = 0.63 eV). The crosses (interpolated by the dotted curve) show the results of 
a calculation using the random phase approximation (RPA). The results of the quantum Monte 
Carlo (QMC) calculations are given with their respective error bars. The region where the system 
undergoes a Mott transition is indicated. The RPA screening remains rather accurate almost all 
the way up to the Mott transition, but fails badly in the insulating region. 



RESULTS AND DISCUSSION 



We have calculated the screening charge An for finite clusters of different size. For 
a given test charge, all calculations give similar results. As an example, Fig. 1 shows 
the screening of a test charge q = 0.5 e as a function of the Hubbard interaction U 
for a cluster of 32 molecules. 

Looking at the results from the quantum Monte Carlo calculations, we can dis- 
tinguish two qualitatively different regions. For weakly interacting systems (U/W 
small) the screening is very good, while it breaks down completely for large U /W, 
where the system is a Mott insulator. RPA somewhat underestimates the screen- 
ing for U/W < 1.0 [7]. For intermediate values of U/W it works surprisingly well 
and actually remains rather accurate almost all the way up to the Mott transition. 
Beyond this point RPA becomes qualitatively wrong. 

The failure of the RPA in the strongly interaction regime can be easily under- 
stood. In RPA it only costs kinetic energy to screen the test charge. Since in 
the strongly interacting region W/U is small, the cost in kinetic energy is small 
compared to the interaction with the test charge. Therefore, when a test charge is 
put on some site, RPA predicts that almost the same amount of electronic charge 
leaves that site in order to minimize the interaction with the test charge. Such 
a calculation neglects, however, the fact that in a strongly correlated system the 
electrons that leave the site with the test charge have to pay a large correlation 
energy. In the QMC calculations this is taken into account properly, leading to the 
observed breakdown in the screening for large U/W. 



The gradual decrease in the efficiency of the screening just before the Mott transi- 
tion that can be seen in Fig. 1 suggests a mechanism for understanding the anoma- 
lous behavior of NH 3 K 3 C6o- At normal pressure this system is an insulator. Under 
pressure it becomes superconducting, with T c increasing with pressure [8]. This is 
in contrast to the other alkali-doped Fullerides, where T c drops with pressure. As 
pressure is applied to the Mott insulator, U/W decreases and the screening be- 
comes more and more efficient, reducing the effective Coulomb interaction /x*, i.e. 
increasing T c . As U/W further decreases, the screening saturates. Now the reduc- 
tion of the electron-phonon interaction (because of the decreasing density of states) 
dominates, and the standard behavior of the alkali-doped Fullerides is recovered. 

To summarize, our results indicate that RPA is surprisingly accurate almost all 
the way up to the Mott transition. Since retardation effects are inefficient, this has 
important implications for the superconductivity in the alkali-doped Fullerides. 
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